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ABSTRACT 



Aims. Molecules that trace the high-density regions of the interstellar medium have been observed in (ultra-) luminous 
(far-)infrared galaxies, in order to initiate multiple-molecule multiple-transition studies to evaluate the physical and chemical 
environment of the nuclear medium and its response to the ongoing nuclear activity. 

Methods. The HCN(l-O), HNC(l-O), HCO + (1-0), CN(l-O) and CN(2-1), CO(2-l), and GS(3-2) transitions were observed 
in sources covering three decades of infrared luminosity including sources with known OH megamaser activity. The data for 
the molecules that trace the high-density regions have been augmented with data available in the literature. 
Results. The integrated emissions of high-density tracer molecules show a strong relation to the far-infrared luminosity. Ratios 
of integrated line luminosities have been used for a first order diagnosis of the integrated molecular environment of the evolving 
nuclear starbursts. Diagnostic diagrams display significant differentiation among the sources that relate to initial conditions 
and the radiative excitation environment. Initial differentiation has been introduced between the FUV radiation field in 
photon-dominated-regions and the X-ray field in X-ray-dominated-regions. The galaxies displaying OH megamaser activity 
have line ratios typical of photon-dominated regions. 

Key words, infrared: galaxies - ISM: molecules - radio line: galaxies - galaxies: active - galaxies: starburst, - Masers 



1. Introduction 

Large amounts of high-density molecular gas play a cru- 
cial role in the physics of (ultra-)luminous (far-) infrared 
galaxies ((U)LIRGs), giving rise to spectacular starbursts 
and possibly providing the fuel for the emergence of ac- 
tive galactic nuclei (AGN). A strong linear relation has 
been found between the far infrared (£fir) and HCN 
luminosity, the latter being an indicator of the high 
density (nfrL)) <T0 4 cm -3 ) component of the interstel- 
lar medium. This arises predominantly from the nuclear 
region and indicates a close relation to the nuclear star- 
burst e nvironment and the pro duction of the FIR lumi- 
nosity (|Gao fc Solomon! l2004al ). The molecular gas con- 
tributes a significant fraction to the dynamical mass of 



clear regions of starburst and AGN nuclei. While much 
of the energy generation of FIR-luminous galaxies has 
traditionally been attributed to AGN activity using their 
optical signatures, the NIR and radi o signatures suggest 
significant starburst c ontributions (jGenzel et al.l 1 1998k 
Baan fc Klocknerll2006f ). Circum-nuclear star formation in 
the densest regions serves as a power source for the ma- 
jority of ULI RGs such that they even mimic the presence 
of an AGN jBaan fc Klocknerl |2006h . The ULIRG pop- 
ulation is a l so characterized by OH megamasers (MM) 
(Baan 119891: iHenkel fc Wilson! fl990h resulting from am- 
plification of radio continuum by FIR-pumped foreground 
gas. Few H 2 CO MM are also found a mong this population 
(jBaan et al]ll993t lAraya et alJl2004h . 



the central regions of FIR galaxies (jHenkel et al 
Young fc Scovillelll99lh . 



1991 



Emissions of molecules that trace the high-density re- 
gions in LIRGs and ULIRGs may serve to study the 
characteristics of the extreme environments in the nu- 
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The observed molecular emissions from the high- 
density components could carry the signature of the nu- 
clear region of the ULIRG and its nuclear power gen- 
eration. Recent studies consider the HCN (1—0) versus 
CO (1—0) relatio n and its implications fo r the star for- 
mation activity ( Gao fc Solomon 2004al lbT). Discussions 
are ongoing in the literature about whether such molcc- 
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ular emissions are accurate tracers of the high-densit 



medium at the site of star formation ( Wu et al.l 1200 



Gracia-Carpio et al~ll2006l : Papadopoulod 



ity 

)5; 



20071) . It should 



state models for higher-density gas at higher tempera- 
tures. An HNC/HCN abundance ratio greater than unity 
suggests a rather quiescent, low temper ature gas. In dark 



bo evident that only multi-transition and multi-molecule 
(multi-dimensional) studies can describe the multi- 
parameter environment of the nuclear activity and further 
establish the connection between the molecular signature 
and the nature of the nuclear energy source. The transla- 
tion or extrapolation from the characteristics of individual 
star formation regions in our Galaxy to integrated regions 
in nearby galaxies at larger distances or even to unresolved 
nuclear regions in distant galaxies requires detailed multi- 
species observations and comparisons combined with elab- 
orate theoretical modelling. The interpretation of molec- 
ular line data has focused on the nature of the local ex- 
citation mechanisms in the form of photon dominated re- 
gions (PDR) with far-ultraviolet radiation fields and X- 
ray dominated regi ons _[XDR) (|Lepp fe Dalgarnol Il996t 
Meijerink fe Spaansl 20051 ). The molecular properties of 
these two different regimes can be used as diagnostic tools 
in interpreting the integrated properties of galaxies. 

This paper reflects a multi-dimensional approach to 
understanding molecular line emissions that started in 
the early 90's, but is only now coming to fruition. This 
study originally aimed to establish a better link between 
molecular megamaser activity and the molecular prop- 
erties of luminous FIR galaxies but it has changed into 
a general study of ULIRGs. Here we present data on 
the molecular characteristics of normal and FIR galax- 
ies across the available range of Lfir- Data of the CO, 
HCN, HNC, HCO+, CN, and CS line emissions ob- 
tained for a representative group of 37 FIR-luminous and 
OH megamaser galaxies has been joined with the ad- 
ditional data of 80 sources presente d in the literature 
Early studies of the properties of CS (IMauersberger et al 
mm, HCO+ (lNguven-Q-Rieu et all Il992h. and HNC 



(jHuttemeister et al.lll995l ) using small numbers of sources 
studied the relations between line and FIR luminosities 
and the nature of the central engine. Our molecular emis- 
sion data show clear tendencies that cover a wide regime 
of nuclear activity. This database forms the basis for a first 
evaluation of the emission line properties and for further 
study and modelling of the line properties. 

2. Characteristics of high-density tracers in 
galaxies 

Extragalactic line ratios do not display the extreme ratios 
in the Galaxy, because they are integrations of ensembles 
of different types of clouds within one beam. Nevertheless, 
clear variations reflect the response of the molecular emis- 
sion to the varying physical and chemical environments 
characterizing the regions of nuclear activity. 

HCN and HNC Molecules - Steady-state excita- 
tion models show that the abundance of HCN and its iso- 
topomer HN C decreases with increasing density and tem- 
perature (|Schilke et al.lll992f ). Values for the HNC/HCN 
abundance ratio below unity are consistent with steady- 



clouds this ratio can be more than 3.0 (Churchwell et al 



19841 ) but in warm environments in giant mole cular clouds. 



the r atio may decrease to values of 0.5 to 0.2 (jlrvine et al.l 
19871 ). In the Orion Hot Core the observed HNC/HCN 
column density ratio can be as low as .013 with val- 
ues of 0.2 on the ridge (jSchilke et al. I ll992l) . The observed 
HNC/HCN intensity ratio is found to vary significantly 
from 0.16 to 2 among galaxies with similar HCN/CO in- 
tensity ratios (lAaltp et al.ll2002l). whi ch may be connected 
with IR pumping (|Aalto et al.l 120071 ) . This suggests that 
HNC is not reliable as a tracer of cold (10 K) high-density 
gas in LIR galaxies, as it is in the Galaxy. HCN has thus 
been used to determine t he gas masses and the related 
star formation efficiencies ( Gao fe SolomorJl2li04bh . 

Our exploratory LVG simulations suggest that, for 
a given density, column density, and temperature, 
HCN(l-O) and HNC(l-O) have similar intensities, but 
that HCO + (l— 0) is weaker. At higher temperatures HNC 
tends to be selectively destroyed in favor of HCN as 
long as the medium is not highly ionized. The canonical 
HNC/HCN abundance ratio is 0.9, if the standard n eu- 
tral production paths are used (jGoldsmith et al.lll98ll ). In 
a highly ionized medium, low HNC /HCN ratios would be 
prevented since H CNH + can form HCN and HNC wit h 
equal probability (|Aalto et alJ 120021 : IWang et al.l l2004h . 
The passage of shocks may selectively destroy HNC and 
could significantl y reduce the stea dy state HNC/HCN 
abundance ratio ( Schilke et al.l 19921). 

The HCO+ Molecule - iHuttemeister et al.l (|l995l ) 
suggest a simplified picture where the molecules trace dif- 
ferent gas components: HCN in the warm and dense part, 
HNC in the cool and slightly less dense part, and HCO + 
in the slightly less dense part, that is penetrated by cos- 
mic rays and/or UV photons. The HCN/HCO + intensity 
ratio may vary from 0.5 to >4, which indicates that the 
mol ecules probe distinct physic al regions in nearby galax- 
ies (|Nguven-Q-Rieu et al.lll992l ) . Alternatively, there have 
to be large-scale abundance differences and/or different 
types of chemistry. The excitation of HCN and HCO + 
in these galaxies suggests low optical depth of less than 
unity and very low excitat ion temperatures of < 10 K 
( Nguven-Q-Rieu et al.lll992h . This conclusion is confirmed 
by studies of many species in N GC 4945 and other nearby 
galaxies fe.g. JWang et al. 20041 ). A multi-transition study 
of the high-density medium in Arp 220 and NGC 6240 sug- 
gests that HCN would serve better as a tracer of dense gas 
than HCO+, which can be more su b-thermally excite d and 
as a radical avoids high densities ( Greve et al. 20061 ). 

Shock propagation can strongly change the gas-phase 
chemistry such that HCO + is enhanced relatively to 
HCN in gas associated with young supernova remnants 
([Dickinson et al.lll980l : IWoottenlll98lh . The shocked region 
of IC 443 shows an enhance ment of two orders of m agni- 
tude in HCO + relative to CO (jPickinson et al.lll980f) . This 
enhancement of HCO + would result from increased ion- 
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ization due to c osmic rays in the shocked dense material 
( Wootten|[l98ll ; lElitzurlll983h . 

The CN Molecule - The abundance of the 
cyanid radical CN is enhanced in the outer regions 
of molecular clouds in PDR dominated environment s 
( Greaves fc Church 1996t Rodriguez-Franco et al. 1998 ). 
At larger depths th e CN abundance declines rapidly 



( Jansen et al. 1995bh . To a significant extent CN is a 



photo-dissociation product of HCN and HNC such that 
e.g. the CN/HCN abundance increases in the outer regions 
of molecular cl ouds that are exposed to ultrav iolet (FUV) 
radiation (see iRodriguez-Franco et al. 19981 ). Therefore, 
the CN/HCN (1—0) line intensity ratios may be used to 
probe the physical and chemical conditions for Galactic 
and extragalactic PDR sources. 

In the Orion bar the CN abundance is enhanced at 
the inner edge of the PD R, but it declines at larger 



depths (jJansen et al.lll995aj ). The Orion-KL region shows 
HCN (1—0) emission that is significantly stronger than the 
HNC and CN (1 —0) emission, which is p redicted for warm 
(T ;> 50 K) gas (|Ungerechts et al.lll997t ). 

Using CN (1-0) and CN(2-1), as well as HCN (1-0) 
measurements across the i nner star-forming molecular 
disk in M82, iFuente et all (|2005h find that Vcn/Vhcn 
w 5. They argue that such a high ratio is indicative of 
a giant and dense PDR bathed in the intense radiation 
field of the starburst and conclude that Ay < 5 in the 
clouds of M82, because for optically thicker clouds the 
CN/HCN column density ratio would be smaller than ob- 
served. The effect of enhanced CN/HCN abundance and 
line temperature ratios is confirmed b y theoretical models 
of CN and HCN gas-phase chemistry (jBoger fc Sternberg] 
20051 ). This suggests that the CN/HCN abundances ratio 
in molecular clouds may be used as probes of FUV and 
cosmic-ray d riven gas-phase chemistry , for a wide range 
of conditions. Boger fc Sternberg! ( 20051 ) find that in dense 
gas, CN molecules are characteristically and preferentially 
produced near the inner edges of the C H zones in the 
PDRs. In addition, chemical models ( Krolik fc Kallmanl 
1983tlLepp fc Dalgarndll996l) indicate an enhancement of 
CN under high X-ray ionization conditions. The promi- 
nence of HNC and CN relative to HCN has been consid- 
ered a me asure of the evolut ionary stage of the respective 
starburst (jAalto et al.l l2002n . This way the objects with 



strong HCN and weak HNC and CN are dominated by a 
warm dense gas environment early in the evolution. 

The CS Molecule - Multi-transition studies of CS in 



selected Galactic clouds by ISnell et al.l (|1984T ) show that 



there is considerable density variation across their core 
regions, because the lower (column) density outer regions 
contribute little to the total CS emission. The cloud core 
is approximately two orders of magnitude denser than 
the extended cloud, and the density rises very steeply 
towards the core. LVG-modelling of the data suggests a 
mean density of 5 to 10 x 10 5 cm~ 3 with cloud core 
masses ranging from 3 x 10 5 to 2 x 10 3 M Q . The CS 
abundance varies considerably, peaking in the clumped 
high-density medium. Chemical studies suggest that the 



fractional abundance of CS i s sensitive to both th e abun- 
dances of sulfur and oxygen ( Graedel et al. 1982j ). 

Although the CS lines may be brightest toward star- 
forming cores in regions of high volume-density, ubiqui- 
tous low-level CS emission, probably due to sub-thermally 
excited low-density gas, is found to dominate th e emis- 
sion from the i nner Galaxy (|McQuinn et al.ll2002f) and its 
central region (jBallv et al.lll988f ). On the other hand, the 
study of the excitation state of CS in the integrated nu- 
clear environment of Arp 220 and NGC 6240 suggests that 
CS probes the densest gas in the high-densit y medium 
better than HCN or HCO+ (|Greve et al.ll2006h . 



3. Observations 

3.1. I RAM Pico Veleta (PV) observations 

Observations with the 30-m iram3 telescope at Pico 
Veleta (PV) (Spain) were carried out in Sept. 1994, May 
1996, Jan. 1997, and May 1997. At the frequencies of the 
A - 3.4 mm J=l-0 transitions of HCN, HNC, and HCO+ 
the full width to half power beam size was 25"; for the A 
= 2.6 and 1.3 mm CO J=l— and 2—1 transitions it was 
21 and 13". Typically, observations were started with the 
CO J=l— and 2—1 lines, then the lower frequency fron- 
tend was tuned to HCN, HNC or HCO + while keeping 
CO(2— 1) as a sensitive indicator of pointing quality. 

During the first session in 1994, the 3 mm SIS re- 
ceiver had a 512 MHz bandwidth and a receiver tempera- 
ture of 80—120 K below 100 GHz (image sideband rejec- 
tion 5-10 dB). Starting in 1996, the renewed SIS receivers 
were single-sideband tuned with image sideband rejec- 
tions of 25— 40 dB. Receiver temperatures were T rec (SSB) 
= 55-85 K at a wavelength of 3.4 mm, 100- 140 K at 
2.6mm, 130K at 2.0mm, and 100-180K at 1.3mm. The 
observations were done under varied atmospheric condi- 
tions, which were best for the runs in 1997. The system 
temperatures during these sessions were 200— 300 K for 
the 80-90 GHz range, 260-400 K around 150 GHz, and 
700-1200 K around 230 GHz. 

During the first session in 1994 we used the filter- 
bank backend (FB) in a 512 x 1MHz mode and the 
acousto-optical spectrometer (AOS) in a 512 x 1 MHz 
mode. During the 1996 and 1997 observations we used 
the AOS in 512 x 2 MHz mode in parallel with the FB in 
256 x 1 MHz mode. The spectral resolution of the observa- 
tions with 1MHz channels was 3.2 km s -1 at 115 GHz and 
1.6 km s" 1 at 230 GHz. 

The observing procedure included the use of the wob- 
bling secondary mirror in order to minimize baseline prob- 
lems, which allowed switching between source and refer- 
ence positions placed symmetrically at offsets of 4' in az- 
imuth. Each wobbler phase had a duration of 2 s. From 
continuum cross scans through nearby pointing sources, 
we estimate the absolute pointing uncertainty to be of 
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Ss3". The positional alignment between the receivers is 
estimated to be better than 2". For all sources single nu- 
clear pointings were used. 

The results of our observations are presented in units 
of ma in beam brightness temperature (T m b; e.g. Downed 
19891 ). Using beam efficiencies of 0.78, 0.75, 0.69, and 0.52, 
we find main beam sensitivities of S/T m b~ 4.86, 4.72, 4.62, 
and 4.44 Jy/K at A - 3.4, 2.6, 2.0, and 1.3 mm. The cali- 
bration is affected by the pointing, baseline stability, and 
receiver stability from session to session. The overall ac- 
curacy of the data is estimated to be 15%. Stability of 
calibration between the four observing periods has been 
ensured by using the same calibration sources during all of 
the sessions. All data has been calibrated using the stan- 
dard chopper wheel method. 



3.2. SEST observations 

Two observing sessions were done at the Swedish/ESCd 
Submillimeter Telescope (SEST) at La Silla (Chile) dur- 
ing January 1995 and January 1997. System performance 
during the first session prevented detection with sufficient 
signal-to-noise ratio for the high-density lines of all our 
target sources except for IR07160— 6215. The data pre- 
sented in this paper have been predominantly obtained 
during the successful second session. 

The observations were done using two Low-Resolution 
AOS systems with 1 GHz bandwidth having a total 
of 1440 channels for a single polarization. The chan- 
nel spacing was 1.4 MHz, corresponding to 5.2 and 3.6 
kins" 1 at 80 and 115 GHz. In 1997 the receiver tempera- 
tures were T rcc (SSB) = 100-180K and 120-140K in the 
80— 116 GHz and 130— 150 GHz frequency ranges, respec- 
tively. Observations were conducted under very favorable 
atmospheric conditions and the actual system tempera- 
tures were between 150 and 200 K in the 80— 90 GHz band 
and between 250 and 350 K around 150 GHz. 

The target transitions for the observations were the 
HCN(l-O) and HCO + (1-0) transitions, which occur m 
a single bandpass, and the HNC(1— 0) transition. Parallel 
observations were done of CS(3— 2) using the other po- 
larization. CO (1—0) observations were done for sources 
without known earlier detections. Single nuclear pointings 
were used for all sources. Pointing observations indicate 
that the absolute pointing uncertainty is of the order of 
&3". 

Calibration was done with a chopper wheel and us- 
ing the T£ intensity scale. The spectral data has been 
converted into a main-beam brightness temperature T m b, 
using a main-beam efficiency of 0.75 at 86 GHz, 0.70 at 
115 GHz and 0.60 at 150 GHz. Conversion factors to flux 
units are 18.75, 18.90, and 20.50 Jy/K. The beam size is 
respectively 57", 45", and 34". 
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4. Results 

The observational results have been presented in Online 
Figures A.l and A. 2 and in Tables B.l and B.2. In follow- 
ing evaluations, the luminosity distance Z?l has been taken 
from the literature for sources with radial velocities below 
about 1500 kms" 1 and determined from known emission 
line velocities or from NED using Hq — 72 kins" 1 and qo 
= 0.5. The FIR luminosity Lfir has been determined from 
the 60 and 100/ zm IRAS fluxes using an I R-color correc- 
tion factor 1.8 ( Sanders fc Mirabell fl99i Solo mon et al.1 
Il997h . 

IRAM Pico Veleta Results - Our Pico Veleta 30-m 
observations of 18 luminous FIR sources have resulted in 
66 detections of emission lines for the following transitions: 
18 for CO(2-l), 11 for HCN(l-O), 8 for HNC(l-O), 11 
for HCO+ (1-0), 8 for CS(3-2), 6 for CN(l-O) and 4 for 
CN(2-1). The target sample for HCN, HNC, and HCO+ 
was selected on the basis of FIR luminosity, the presence 
of OH MM activity, or the strength of the CO J=l-0 
line obtained earlier with the 12-m NRAO telescope (Baan 
& Freund 1992, unpublished). Sources for observations of 
the cyano radical CN were selected on the basis of the 
strength of the CO and HCN lines. The CS(3-2) data 
has been obtained in parallel to the observations at lower 
frequencies. High-density tracer transitions in sources with 
published detections at the time of the observations were 
not re-observed. Individual spectra toward the nuclear po- 
sitions obtained during our observations are presented in 
T m b units in Figs. A.l and A. 2. Results from gaussian fit- 
ting are presented in Table B.l after further conversion 
into units of Jy kms — . 

SEST results - The observations of 20 sources with 
the SEST 15-m telescope have resulted in a total of 52 new 
detections of emission lines in 16 sources for the follow- 
ing transitions: 12 for CO(1-0), 11 for HCN(l-O), 8 for 
HNC(l-O), 13 for HCO+ (1-0), 7 for CS(3-2), and 1 for 
CN(1— 0). Target sources during these sessions have been 
select ed on the basis of the strength of known CO emission 
lines ([Mirabel et alJll99(A iGaray et al.lll993h . and of FIR 
luminosity taken from the IRAS catalogue. The source 
NGC 660 is common between the Pico Veleta and SEST 
data sets. High density tracer lines in sources with pub- 
lished detections at the time of the observations were not 
re-observed. Individual spectra toward the nuclear posi- 
tions are presented in Fig. A. 3 in T m b units and the results 
from gaussian fits arc given in Table B.l after conversion 
into units of Jy kms" 1 . Target sources without any line 
detection have been omitted from the tables. 

Detections from the Literature - The detections 
of our combined sample has been augmented with detec- 
tions of molecular lines tracing the high-density compo- 
nent reported in the literature. In addition, we have used 
these detections to verify and cross-calibrate our own de- 
tections where possible. Thus we obtained detections of at 
least one of the HCN(l-O), HNC(l-O), HCO+(1-0), CN 
iV=l-0 and N=2-l, and CS(3-2) transitions for a total 
of 80 additional sources. A total of 167 individual pub- 
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lished records have been used, w hich includes a significan t 
number of HCN detections from lGao fe Solomon! (|2004a! ). 

The Integrated Source Sample - Our own sam- 
ple of 37 sources has been augmented with the data of 80 
sources from the literature. Table B.l presents the avail- 
able data on molecules that trace low-density and high- 
density components of the medium. The database con- 
tains a total number of 202 reliable detections of high- 
density tracers, and 110 for CO (1—0) with 32 accompa- 
nying CO (2—1) line detections. This (still) represents an 
incomplete record for many of the sources. The database 
has 84 detections plus 18 upper limits for HCN, 28 + 
13 for HNC, 42 + 10 for HCO+, a total of 19 + 6 for 
CN(l-O), 13 + 7 for CN(2-1), and 17 + 20 for CS(3-2). 
For 23 sources a complete set of HCN (1-0), HNC (1-0) 
and HCO + (l— 0) tracer transitions were obtained and 10 
sources have one out of three as an upper limit. 

All data records for our integrated sample have been 
converted from T m b into units of Jy kms -1 for the PV 
and SEST results or from antenna temperature for other 
telescopes. We also note that some data in the literature 
has been obtained with various single-dish instruments: 
the 30-m Pico Veleta telescope, the 15-m SEST telescope, 
the 12-m NRAO telescope at Kitt Peak, the Onsala 20- 
m telescope, and the 14-m FCRAO telescope. The inclu- 
sion and comparison of this data has been complicated 
by the incomplete reporting of observing and calibration 
procedures. In many cases the CO (1—0) and HCN(1— 0) 
line strengths from different data records have been used. 
General agreement was found for line data from multiple 
source records with errors within 30% . Some publishe d 



line data records, particularly in lGao fc Solomon! (|2004a! ), 
were not compatible with our data and other available 
data and could not be used in the compilation. As a re- 
sult some line ratios are also different from those in other 
publications because of inconsistencies between our own 
data and some published data. The source NGC 660 was 
included in both our observing samples and a comparison 
of the HCN and HCO + lines shows agreement within 10%. 
We also note that some of our own lower quality spectra 
and their unreliable fits have not been used in the evalua- 
tion and have been designated as upper limits in our data 
base in Tables 1 and 2. 

Evaluating the Data - The collected data has been 
obtained at different observing beams and different beam 
filling factors. For a given telescope, the beam sizes at 
the observing frequencies of HCN, HNC, and HCO + are 
very similar and they differ by some 25% from those of 
CN(l-O) and CO (1-0). Since the emission of the high- 
density tracers is well-confined to the nuclear regions of 
the galaxies, it will generally not fill the beam. On the 
other hand, the CO (1—0) emission and possibly also the 
CO (2—1) emission will be more extended and could fill 
the beam for nearby galaxies. The extent of a typical 
CO (1 —0) emission region i n ULIRGs varies between 2 and 
3 kpc (jDownes fc Solomonlll998f ); this would result in an 
underestimate of the CO emission at radial velocities less 



than 2000 kms for the PV telescope and less than 1000 
kms -1 for the SEST telescope. 

Because of their similar frequencies, the line ratios of 
the high-density tracers are not affected by different beam 
filling factors. Their ratios with CO (1-0) would only be af- 
fected by an underestimate of the total CO emission. The 
line intensity ratio between CO (1-0) and the high den- 
sity tracers used throughout this paper will therefore be 
representative for more distant galaxies, while they could 
be lower limits for nearby galaxies. We will return to the 
effect of relative filling factors in the evaluation of the CO 
line ratios below (Sec. [57 



5. The CO line characteristics 

5.1. Line luminosities 

The J=l— and 2—1 mm-wave lines of CO observed by 
us are easily saturated and thermalized, such that they 
are likely not tracing the high density regions (n(H2) > 
10 4 cm -3 ) of the ISM where the star formation occurs. The 
CO characteristics of our sample sources are shown in Fig. 
[TJ Slopes (indicated in the figures with s) have been fitted 
to t he data using a standard Le ast Squares Fit method 
(see Bevington fc Robinson 19921 ). To evaluate the signifi- 
cance of the fits, the correlation coefficient (indicated by r) 
has also been calculated. For these calculations upper lim- 
its have been omitted, and all data have been given the 
same statistical weight, since only few data points have 
known errors. Both CO transitions display a clear depen- 
dence on the FIR luminosity with slopes of 0.74 ±0.04 
(110 sources; r = 0.87) and 0.91 ± 0.10 (32 sources; r = 
0.85). At high luminosities the scatter in Leo is a factor 4 
but increases at lower luminosities to a factor of 10. This 
scatter can be attributed to variations of the star forma- 
tion efficiency (SFE) Lfir/M(R2), ranging between 1 and 
100 L0/M0, and (in part) to underestimates of the CO 
luminosity for the lower luminosity (and nearby) sources, 
where single-dish instruments only sample the central re- 
gion. 

The slopes for both CO curves are significantly less 
steep than unity, possibly related to optical depth ef- 
fects or to an increase in the star formation efficiency 
SFE at high Lfir- Assuming constant Loo-to-M (H2) con- 
versions, earlier studies suggested a linear Lfir - A/(H 2 ) 
relation with a slope that steep ens at high Lftr due t o 
a strong variation of the SFE (|Young fc Scoville![T99l . 
Similarly iBraine fc Combesl ( 1992^ present a linear rela- 
tion between the CO(2-l) and FIR data for L F ir < 10 10 5 
L/0 . The relation with slope less than unity deduced from 
our data for Lfir > 10 10 L Q is consistent with character- 
istics of the upper luminosity range of the earlier data. 
Our fit would suggest SFE oc Lfir ' 25 for the higher end 
of the Lfir distribution, or an equivalent variation in the 
Lco-to-M (H2) conversion factor. 
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Fig. 1. Integrated line luminosities versus FIR luminos- 
ity for the integrated sample of sources, a) CO (1—0) lu- 
minosity versus FIR luminosity, b) CO (2—1) luminosity 
versus FIR luminosity, and c) CO (2-l)/CO (1-0) lumi- 
nosity ratio versus FIR luminosity. There are 110 reliable 
data points for CO(l-O) and 32 for CO(2-l). The fitted 
slope (s) and the correlation coefficient (r) are given in the 
lower right corner of each frame and the fit does not take 
into account any displayed upper limits. The data points 
in the diagrams are either squares for reliable values or 
triangles for upper and lower limits. Filled symbols indi- 
cate the 24 OH or H2CO MM sources in the sample. Open 
symbols may also represent undetected OH MM because 
of beaming effects or weak maser action with similar ISM 
conditions. 



5.2. The CO line ratio 

The 7co(2— l)//co(l — 0) line ratio has been displayed 
for 32 sources of which five have values greater than 
2 (OH MM IRAS 15107+0724, NGC660, NGC3256, 
IRAS 06259-4708, and IRAS 18293-34 13; Fig. EDO - 



Earli e r results from nearby s piral galaxies (jRadford et al 



199lt iBraine fc Combes! 11992 ) as well as the data from 
this study show a concentration of data points up to 1.5 
with a few outliers up to three. 

An observed line ratio of > 1 for a homogeneous gas 
could imply that the gas is optically thin, which is not 
realistic for the nuclear environments considered in our 
sample. An optically thick cloud model with external heat- 
ing of small dense clouds, i.e. with hot young stars blow- 
ing holes into the molecular structures, can make the 
CO (2—1) line stronger than the CO (1—0) line, if the tem- 
perature decreases by about 30— 60 K (factor 2) from the 
r = 1 surface for the CO (2—1) line to the r = 1 sur- 
face for the CO (1-0) line (jBraine fc Combeslll992r i. The 
higher opacity of the CO (2—1) line as compared with the 
CO (1—0) line, in almost all scenarios, leads to a higher 
beam filling factor for the CO (2—1) line and thus a line 
intensity ratio larger than unity. 

Early evaluation of the Ico to hydrogen mass MjB.^) 
conversion factor in galaxies by iMalonev & Black! (|l988l ) 
assumed thermalized CO rotational levels with T cx = Tkin- 
This would predict CO (2-l)/CO (1-0) brightness tem- 
perature ratios in the range of 0.9 to 1.1 in the Galaxy. In 
order to explain similar extragalactic data. iRadford et al 



([19911) show that sub-thermal excitation (T cx < Tki n ) in 
the n(H2) = 10 2 to 10 3 cm -3 range could result in inte- 
grated CO (2-1) /CO (1-0) line ratios in the 1.4-2 range, 
and beam-corrected (intrinsic) 7t>(2— l)/2~b(l— 0) ratios in 
the range of 0.6—0.75. 

An evaluation of the CO properties requires considera- 
tion of the relative filling factor of the two CO transitions 
(Sec. SI). Here the conversion from the observed T m b ratio 
to an intrinsic Tb ratio involves an antenna beam cor- 
rection factor (g? 2 + &2i)/(^s + &io)> which expresses the 
relative size of the source d s and the beam size b at the 
two transitions. For a representative size of 2 kpc for the 
prominent CO emission regions in luminous FIR galax- 
ies, the intrinsic Tb ratios will be reduced by a factor of 
1.15 - 2.46 (for velocities of 750 to 7500 kms" 1 ) relative 
to the corresponding T m b ratios. This range of values sug- 
gests that most of the observed ratios could have intrinsic 
Tb ratios around unity, which is consistent with a sub- 
thermal excitation o f the low density molecular medium 
(jRadford et al.lll99l! ). Our most extreme value of 3.0 for 
the OH MM source IR 15107+0724 at 46 Mpc would give 
a brightness temperature ratio of 1.4. 

The conversion of integrated CO (1—0) lines to H2 
column densities is expressed as -/V(H 2 )/ico(i-o) = a 
([Strong et al.l ll988). These authors suggested a value a 
= 2.3 1 20 molecules cm _2 /(K kms -1 ), using the proper- 
ties of self-gravitating virialized molecular clouds of high 
optical depths. However, this value is not uniformly appli- 
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cable for the Galaxy and for extragalactic nuclei. Studies 



of th e Galactic Center using COBE data (Sodroski ct al 
1995) and C 18 and the mm dust emission (jDahmen et al 



19981) find that the conversion factor is a highly vari- 



able value, resulting in an overestimate of the mass in 
the nuclear region and an underestimate in the outer 
disk. Studies of the central regions of th e nearby galax- 
ies NGC253 (IMauersberger et al.1 fl996bh and NGC4945 
( Mauersberger et al.l 1996a ) show an overestimate by a 
factor of 6—10. Studies of ULIRGs show that most of 
the CO luminosity comes from relatively low-density gas 
that is d ynamically significant and is located in a molecu - 
lar disk ( Solomon et al.lll997t iDownes fe Solomon! Il998l ) . 
A conservative estimate of the CO— to— H2 conversion 
M(R 2 )/L' co = 0.8 M /(K km s" 1 pc 2 ) for such extra- 
galactic nuclei is based on an a factor about 5 times lower 
than the value for self-gravitating clouds. This conversion 
ratio may be used to convert the CO (1—0) values listed 
in Table B.l into H2 masses. 



6. High-density tracer molecules 

6.1. Line luminosities 

The diagrams of line luminosity against the infrared lumi- 
nosity for all molecular transitions considered (Figs. Q] HI 
and [3]) show a well-behaved distribution with a consider- 
able width. The spread of the data points may be caused 
by a variety of effects, such as non-radiative heating of 
the gas, differences in the initial abundances, and the evo- 
lutionary changes of the ISM resulting from a starburst 
and the depletion of the high-density material (see Sec. [7] 
below). Data points at low luminosity that lie well below 
the center of the line luminosity distribution may result 
from using only the central pointing for nearby galaxies. 

The HCN, HNC and HCO+ Molecules - The line 
luminosities of the source sample of the prominent tracer 
molecules are presented against the FIR luminosity in Fig. 
[2 The slopes for these three molecules are close to unity 
such that HCN (84 sources) has s = 0.98 ± 0.06 (r = 
0.89), HNC (28 sources) has s = 1.18 ± 0.13 (r = 0.88), 
and HCO + (42 sources) has s = 1.10 ± 0.06 (r = 0.95). 
The tracer luminosities provide a slightly tighter relation 
with Lfir than the CO(l— 0) data as displayed in Fig. [T] 
Fits with slope 0.9 ± 0.1 for HCN versus HCO+ and with 
slop e 1.4 ± 0.4 for HCN versus CO(l-O) have been noted 



slop e 1.4 ± u.4 tor iii;rN yersus u uii-U) nave been noted 
bv lNguven-Q-Rieu et all (|l992f ). A slope of 0.97 ± 0.05 



was found for HCN versus I/fir (|Gao fc Solomon! l2004af) 
and slo pes of 0.88 ± 0.15 for HN C versus both HCN and 
HCO + (|Huttemeister et~aT] Il995h . The scatter of the dis- 



tributions is about one decade in the main body of data 
points. 

The CN Molecule - The mm— wave spectra of 
the cyano radical CN have more complex characteris- 
tics than those of the other observed molecular species. 
The CN A^=l— line observations at Pico Veleta have 
been taken at the strongest feature, the J=3/2— 1/2 tran- 
sition at 113.491 GHz that contains one third of the 



total intensity in the optically thin limit under Local 
Thermodynamic Equilibrium (LTE) conditions. For the 
SEST data for NGC 1068, besides the J=3/2-l/2 feature 
a second J— 1/2— 1/2 feature is seen at 113.191 GHz at a 
relative strength of 0.45 (Figs. A. 2). This relative strength 
is close to the LTE value, which could sugges t that the 
lines are optically thin (see IWang et alJliool l. The CN 
A^=2— 1 transition data focus on the J=5/2— 3/2 line at 
226.659 GHz that represents one sixth of the total inten- 
sity (r <1; LTE). The CN dat a has been augmen ted with 
data for both transitions from Aalt o et al. I (|2002h . 

Most of the CN detections are single broad features 
centered at the systemic velocity of the galaxies (Figs. 
A.l-2). Only Arp 220 and IR05414 show double N = 
1—0 and N = 2—1 CN components that cannot be at- 
tributed to other spingroups but only to the two nuclei of 
their respective systems. 

The CN line luminosities in Figure [3Kb show a distri- 
bution that is similar to that of other tracer molecules. A 
slope close to unity is found for CN(1— 0) (19 sources): s 
= 1.03 ± 0.13 (r = 0.89) and for (CN(2-1) (13 sources): s 
= 1.14 ± 0.17 (r = 0.89). The scatter of the distributions 
is less than one order of magnitude. 

The CS Molecule - Our CS (3-2) data from the PV 
30-m and SEST telescopes ( Figs. A.l— 2) have be e n aug - 



mented with detections f rom | Mauersberger et al 
IMauersberger fc Henkell ( 19891 ) . and IWang et al. 



1989), 
2004). 



An early assessment of the CS emission in nearby galax- 
ies shows a clear re lation with Lqq and Lfir an d the 
star formation rate ( Mauersberger fc Henkell 1 19891) . The 
CS line luminosity (17 sources) versus the FIR luminos- 
ity has a higher slope, s = 1.56 ± 0.14 (r = 0.95), as 
comp ared with a slope of unity f or the much smaller sam- 
ple of Mauersberger et al. ( 1989h . This larger slope of CS 
as compared to the other molecules under consideration 
indicates a stronger environmental dependence (Fig.[3fc). 



6.2. Evaluating the relations 

The well-known Kennicutt-Schmidt laws relate the 
surface- density of star formation (Esfr) and the surface- 
density of the gas (Y, gas ). On the basis of HI, CO, and 
Ha data for a l arge sample of spiral and starburst galax- 
ies, [Kennicutt ( 19981 ) finds that the Esfr — E gas relation 
has an exponent a=lA (Esfr cx ^g as )- Analogously the 
relation of the integrated luminosity of the high-density 
gas component involved in the star formation process 
and the integrated star formation rate (SFR) expressed 
as Lftr could also behave as a Kennicutt-Schmidt law. 
Krumholz fc Thompson! (|2007l ) made model predictions 
for the Kennicutt-Schmidt laws assuming a constant star 
formation efficiency. They find that the line luminosity for 
molecules with a critical density lower than the average 
density of the system rises linearly with density. For lines 
with critical densities larger than the median gas density 
the luminosity per unit volume increases faster than lin- 
early with the density. The star formation rate also rises 
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super-linearly with the gas density, and the combination 
of these two effects produces a close to linear correlation 
(a ~ 1) for the high critical density molecules and a super- 
linear relation for the low density tracers (a ~ 1.5; the 
classic Kennicutt-Schmidt law). They also predict that 
generally there may be more intrinsic scatter in this rela- 
tion for low critical density molecules than for high density 
tracers. 

Our data mostly confirms the findings of 
Krumholz & Thompson! (|2007l ). Our CO distribu- 
tions (Fig. [1} show indices of a=1.40 and 1.23 (the 
inverse of the slope s in our diagrams) and have a 
relatively large spread. The relations between Lfir, and 
the line luminosities of high density tracers HCN(1— 0), 
HNC(l-O), HCO + (1-0), CN(l-O), and CN(2-1) are 
all linear within the errors (see Figs. [5] and [3]). However, 
the relation for CS (3—2) is sub-linear with a=0.66. Since 
CS has a critical density close to that of HCN and higher 
than that for CO, CS is expect ed to have a linear relation 
(jKrumholz fc Thompsor] 120071 ). This may indicate that 
this species is influenced by additional chemical and 
environmental effects. 



6.3. Tracer intensities relative to CO (1—0) 

In Fig. 0] and \5\ the integrated intensities for the to- 
tal sample are compared with the CO(l— 0) emission 
that represents the larger scale low-density component. 
Ignoring NGC7469, which seems to have unusual high 
ratios, and upper limits, the emission line ratio varies 
from 0.01 to 0.18 for HCN, from 0.01 to 0.08 for HNC, 
from 0.01 to 0.21 for HCO+, from 0.01 to 0.15 for 
CN, and from 0.01 to 0.11 for CS. An early study of 
CS suggests a relative line rat io in the 0.01—0.08 range 
(Mauersberger fc Henkellflfjfl . The upper range for the 
HCN/CO intensity ratio o f our current sample is com- 
parable with the sample of Gao fc Solomonl ( 2004al) . The 
range for HNC and HCO + is found to be (slightly) smaller 
than for HCN (when ignoring the outliers/upper limits) 
(see also lGracia-Carpio et al.ll2006l ). 

We note that the distribution of characteristic line ra- 
tios for all molecules increases with FIR luminosity, which 
gives an upwardly curved lower boundary for the distri- 
bution at higher Lfir- The highest values for all ratios 
are found at Lfir > 10 10 5 L Q . As discussed in Sec. [7] 
below, a decreasing line strength of a high-density trac- 
ers as compared with CO during the course of a nuclear 
starburst would be the signature of the depletion of the 
nuclear high-density component. 



7. Modelling of a ULIRG outburst 

7.1. The outburst model 

The FIR light curve of an ULIRG is determined by the ra- 
diative energy generated by a (circum-)nuclear starburst 
possible in hybrid mode with an AGN, that is re-radiated 
by a torus-like dusty environment. During such an FIR 



outburst, the low-density and high-density molecular com- 
ponents of the nuclear region evolve differently. The large- 
scale signature of the low-density gas, as represented by 
CO (1—0) emission, would only be partially altered by the 
centralized star formation process. On the other hand, the 
centrally concentrated high-density gas forms the struc- 
tural basis for the star formation process and will be slowly 
consumed, destroyed or dispersed by the newly formed 
stars. 

In a simplified scenario, the instantaneous FIR lu- 
minosity of the ULIRG represents the rate at which 
the nuclear activity generates FIR energy and con- 
sumes/destroys its high-density molecular material. In the 
absence of a representative FIR light curve with different 
timescales for the rise and fall of the FIR luminosity, we 
choose a diffusion-like expression as a response to a sudden 
start of the star formation. This function may be defined 
for t > as: 



Lfir(£) = 1-35 L Fm (max) 



2.-> 



-T/t 



(1) 



where LFiR(max) is the maximum luminosity and T is 
the characteristic timescale of the burst. The duration of 
the burst (down to 1% of the peak) is on the order of 10 
x T. A diffusion curve may not be the most appropri- 
ate representation of a FIR light curve but it resembles 
those obtained with starburst-driven models of the FIR 



evolu tion and color properties of ULIRGs (jLoenen et al 
l200fih . 

The initial line strength of the high-density HD(t = 0) 
component can be defined as a fraction /3o of that of LD, 
the larger-scale low-density component, that is assumed 
to remain unchanged during the outburst. As a result, the 
line ratio of the high- and low-density components varies 
with time during the FIR outburst as: 



HD(t) 
LD 



00 



i-r 



/„' L Fm (t)dt 



(2) 



where T is the total fraction of the initial (3$ of the HD 
emission component that is destroyed (disappears) dur- 
ing the whole course of the outburst. The instantaneous 
HD component depends on the fraction of the FIR energy 
generated up till time t and the total FIR energy radiated 
during the outburst LpiR^nt- 

7.2. Modelling results 

The results obtained in these simulations are presented in 
Fig. [6] for three FIR light curves with peak luminosities 
log LYYn{max) = 11, 12, and 12.5 in order to cover the 
range of observed FIR luminosities. The FIR light curves 
(upper frame) are used to determine the evolving lumi- 
nosity of a high-density component (middle frame), and 
the evolving emission ratio of the high- and low-density 
components (bottom frame). A light curve having a more 
delayed peak would give more rounded curves and the 
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turn-around point would occur at higher HD(t)/ LD val- 
ues. 

This simple time-dependent scenario already shows 
some distinctive behavior of the high-density component 
during a FIR outburst, that is relevant to the high-density 
tracer data presented in this paper: 

1) The rapid rise segment (indicated as (a) in Fig. [6]) of 
the star formation process and the resulting FIR outburst 
ensure that the data points show rapid evolution and are 
relatively sparse at the onset of the curves. During the 
peak of the curve the depletion rate of the high-density 
component is highest. After the luminosity peak the data 
points in the decay segment of the three curves (indicated 
as (b) in Fig. O will be closer together. Because of the 
increased probability of finding a source in this state, our 
source samples will be dominated by sources in the inter- 
mediate and later stages of evolution. 

2) The luminosity of the high-density component versus 
Lfir will display a separation between the upward (a) and 
downward (b) luminosity segment due to the ongoing de- 
pletion (middle frame Fig. [6]). Sources with the highest 
peak luminosities and the largest depletion factors follow 
the outermost evolutionary track and will end on the lower 
boundary of the distribution. This effect could add signif- 
icantly to the scatter in the diagrams of the high density 
tracers versus Lfir- 

3) The high- to low-density emission ratio versus FIR 
(bottom frame Fig. [6|) reflects the phases of the FIR light 
curve: a rapid horizontal track for the rapid rise segment 
(a) of the FIR luminosity with few data points, followed by 
a downward curved track with accumulation of data points 
at the lower boundary for the decay segments (b). After 
finishing the FIR light curve, a third buildup phase may 
take place with a steady re-building of the high-density 
molecular component in the central regions. A molecular 
buildup would prepare the galaxy for the next star forma- 
tion outburst. 



7.3. Comparison with observations 

The observed relations between molecular line luminosi- 
ties and FIR luminosities in Figs [1] [2] and [3] are in agree- 
ment with the simulations in the middle frame of Fig. [6] 
The observed scatter in the diagrams would partly result 
from differences in the rise and decay stages of an out- 
burst. Furthermore, the track of the highest luminosity 
sources and OH MM would lie at the outer edges of the 
distribution. Similarly, the distribution of data points in 
an Lhd - Lyir diagram could display a steepening to- 
wards higher Lfir values (see Fig. [2J. 

The observed intensity ratios of tracer molecules rel- 
ative to CO (1—0) show organized distributions (Figs. [4] 
and [5j, that are in agreement with the simulations in 
Fig. [6b. The upper regions of the distributions are indeed 
sparsely populated as they represent (near-)initial values 
for sources. From the diagrams it follows that the initial 
fraction (3 of the high-density component would be in the 



0.08—0.18 range, and the depletion rate T for our sample 
of sources would lie in the 0.6—0.9 range. Further evolu- 
tion of the outburst during the decay phase would result 
in the characteristically curved lower boundary observed 
in the diagrams. The highest FIR luminosity sources that 
have passed beyond their peak are also found at this lower 
boundary. 

The distribution of OH MM (and H 2 CO MM) are quite 
conspicuous in the diagrams. Some occupy the sparse re- 
gion of rapid evolution at the highest HD/LD intensity 
ratios (high /3 ), while others have already evolved to the 
curved lower edge of the distribution, as seen clearly for 
the HCN/CO ratio in Fig. Hi. 

8. Molecular excitation environments 

The emission properties of Galactic and also extragalac- 
tic emission regions have been studied and interpreted by 
using two distinct excitation regimes: X-ray Dominated 
Regions (XDRs) and Photon Dominated Regions (PDRs) 
with far- ultraviolet radiation fields. These environments 
are quite distinguishable using the emission characteristics 
of high-density tracers in the cloud cor e regions and low- 
density tracers at the cloud surfaces (iLepp k. Dalgarnol 



ll996l : lMeijerink k Spaansl2005l : lBoger fc Sternbergll2005h . 
The abundance and emission characteristics of the molec- 
ular and atomic lines depend strongly on the available (at- 
tainable) column density submitted to the radiation field 
and the strength of the radiation field. 

Global environmen tal conditions - The PDR 
and XDR mod e lling of Meiierink &: Spaand ( 2005 ) and 
Meiierink et al ] (120061 . l2007h incorporate the physics and 



the chemistry in the local environment using many 
molecules and chemical reactions under local thermal bal- 
ance. A large grid of density and radiative conditions was 
used to predict the cumulative column densities taking 
into account the radiation transfer in a representative slab- 
model incorporating XDR and PDR environments. This 
grid predicts column densities in all relevant energy states 
for the dominant molecules. All high density tracer emis- 
sions are predominantly produced by regions with high 
hydrogen column densities. Relying on the predicted col- 
umn densities for such a gas component, we restrict our- 
selves to a first order analysis of the observed line ratios 
using the ratios of calculated column densities. 

The characteristics of both the PDR and XDR models 
depend on the parameterization of the radiative energy 
deposition rate per hydrogen nucleus P — Ffuv/^h or 
Fx/riu and suggest the following regimes: 

1) For PDRs with high values of P > 5 x 10 -4 erg cm 
s _1 the column density ratios vary for all column densities 
A H > 10 21 ' 5 cm" 2 as: 

A(HCO+)/A(HCN) < A(HNC)/A(HCN) < 1 

while for smaller values for Ah, the ratios switch around 
such that: 

A(HNC)/A(HCN) << 1 < 7V(HCO + )/7V(HCN) 
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and for lower P : 
ratios vary as: 



x 10 6 erg cm s 1 the column density 



1 < N(HNC)/jV(HCN) w iV(HCO + )/7V(HCN) 

2) For XDRs the radiation penetrates deeper into the 
cloud volume than for PDRs. For stronger XDR radiation 
fields (P > 5 x 10 -4 erg cm s _1 ) and for column den- 
sities above the transition point Nh w 10 22 5 cm -2 , the 
observable molecular column densities have ratios: 

1 < 7V(HNC)/7V(HCN) < iV(HCO + )/A(HCN), 

and below this A^h transition point: 

iV(HNC)/iV(HCN) << iV(HCO+)/iV(HCN) < 1, 

while for lower P « 5 x 10~ 6 erg cm s _1 the column 
density ratios vary as: 

7V(HCO+)/iV(HCN) « iV(HNC)/7V(HCN) w 1 

3) For PDR conditions the column density ratio 
AT(CN) /AT(HCN) varies from 0.5 to 2.0 over the respective 
density range of 10 6 to 10 4 cm~ 3 . 

4) For XDR conditions the column density ratio is 
iV(CN)/7V(HCN) >> 1 and may become as large as 10. 

5) For PDR conditions at N H > 10 22 cm" 2 the col- 
umn density ratio iV(CS)/iV(HCN) > 1. At lower column 
densities the ratio may become < 1, which makes CS a 
column density indicator. 

6) For XDR conditions for the whole column density 
range at low values of P ~ 5 x 10 -6 erg cm s" 1 the column 
density ratio is 7V(CS)/jV(HCN) > L At higher values of 
P the ratio remains > 1 for TVh > 10 23 cm~ 2 and is < 1 
for lower iVjj, again making it column density sensitive. 

Evolutionary environmental conditions - The 
characteristics described above are representative of a 
steady state nuclear excitation environment. Realistically 
the behavior of molecules tracing the high-density compo- 
nent is affected by changing evolutionary and environmen- 
tal effects that influence their abundance and the exciting 
radiation fields. 

Assuming that FIR luminous sources are mostly pow- 
ered by (circum-) nuclear starbursts, the evolving star for- 
mation will generate supernovae with shocks enhancing 
the cosmic ray flux, and produce new generations of X- 
ray binaries. Because these star formation products alter 
the excitation environment, an environment initially dom- 
inated by PDR conditions could change into one with a 
combination of PDR and XDR condi tions during the later 



stage s of the outburst (for M82 see: iGarcia-Burillo et al 



200a ISpaans k. Meiierinkj 120071 ). This will modify the 
line ratios from their steady state conditions discussed 
above. For instance, increased cosmic ray production in 
shocks would enhance the relative HCO + abundance while 
shocks also would decrease the relative HNC a bundance 
(|Woottenlll98~0 : lElitzur|[l983 ISchilke et~aDll992h . We will 
refer to this evolved starburst stage as being in a "late- 
stage" as compared with an " early-stage" starburst where 
these effect do not yet play a role. 



During the late-stage starburst also the density and 
temperature structure of the ISM and its clumping would 
be affected, resulting in chemical changes such as the selec- 
tive destruction of HNC in favor of HCN at higher tem- 
peratures (see Sec. [2]). In addition, a certain fraction of 
the high-density medium will be used/destroyed by the 
star formation process. As a result the line emission may 
originate in regions of decreasing column densities, which 
modifies the observed line ratios. 



9. Diagnostics with tracer lines 

In the absence of multi-transition studies using LVG sim- 
ulations involving higher J transitions, we make the as- 
sumption that the integrated line ratios have a diagnostic 
value that is equivalent to the ratios of column densities 
for the observable molecular species. The global charac- 
teristics of molecular column densities described in the 
previous section for XDR- and PDR-dominance may thus 
be used for a first diagnosis of the excitation environments. 

9.1. High-density tracer ratios versus Lfir 

Diagrams of relative intensities for HCN, HNC, HCO + , 
CN, and CS are presented in Figs. [7] and [8l Horizontal 
lines for a ratio equal to unity (and two for CN) have been 
added to the diagrams in order to test a first order divi- 
sion between PDR and XDR line ratios as derived in the 
previous section. However, we find that the unity value is 
not a reliable boundary for ratios with HCO + . For CS this 
line separates high and low column density. The line ratios 
are predominantly determined by regions having the h igh- 
est molecular column densities ([Meiierink et al.ll2007f ) . All 
diagrams show some systematic behavior in that the OH 
MM sources cluster together and that there could be a 
wedge-like upper boundary for the sources at high lumi- 
nosities. 

HCN, HNC and HCO+ lines - The observed ratios 
of HCN, HNC, and HCO+ in Fig.[7]show a range of values 
between 0.2 and 2.0. Modelling results indicate that the 
HCO + /HCN ratio in most PDR environments is expected 
to be smaller than unity, while in XDRs the HCO+/HCN 
line ratio becomes > 1 for the higher column densities. 
The PDR sources would lie in the lower parts (< 1) of the 
upper two diagrams and in the upper part (> 1) of the 
bottom diagram of Figure [7l 

The observed HNC/HCN data shows significant clus- 
tering at 0.5 across the whole Z-fir range, which could 
indicate that the HNC/HCN ratios have been lowered 
by the evolving starburst excitation. In the HCO + /HCN 
diagram in Fig. [7Jd, the OH MM sources are well con- 
fined to the 0.3 and 0.8 range with only IC 694 (Arp299A, 
IR 11257+5850A) as the exception. Rather than finding a 
steady dec rease of the observed HCO + /HCN with increas- 
ing Lfir (|Gracia-Carpio et"aL 20061 ). our larger sample 
presents a more complicated picture that hints at (sys- 
tematic) evolutionary influences on the nuclear medium. 
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Fig. 9. The integrated line ratios of HCN(l-O), HNC(l-O), and HCO + (1-0) versus each other, a) Integrated 
HCO+ (l-0)/HCN(l-0) versus HNC (1-0)/HCO + (1-0) ratios, b) Integrated HCO+ (l-0)/HCN (1-0) versus 
HNC (l-0)/HCN(l-0) ratios, and c) Integrated HNC(l-0)/HCN (1-0) versus HNC (1-0)/HCO + (1-0) ratios. The 
P and X and the dotted lines mark the regions of the line ratios of predominantly photon-dominated regions (PDRs) 
and predominantly X-ray-dominated regions (XDRs) using only the HNC/HCN ratio. The symbols are described with 
Figure 1. 



The diagrams of Figure [7] display systematic behavior 
related to the HCO + and HNC characteristics. The group 
of megamaser sources (mostly) lies to the right of the ver- 
tical line of Lfir > 10 11 L0 but they are in the PDR 
regions of all three diagrams. They are mostly high Nn 
and pure PDR sources. On the other hand, the group of 
sources left of that vertical line lies in the PDR regime for 
the HNC/HCN ratio, straddles the XDR-PDR dividing 
line for the HCO + /HCN ratio, and lies in the XDR region 
for the HNC/HCO + ratio. This effect can be attributed 
to HCO + enhancement and HNC depletion during a late- 
stage starburst phase, that pushes PDR sources to a first 
order XDR regime in the two bottom diagrams. This late- 
PDR behavior is exemplified by the late starburst M82. 

The lowest values of HNC/HCO+ for high 
Lfir sources are for OH MM IC 694 and the two 
AGN sources NGC 1068 (H 2 MM) and NGC6240. 
NGC6240 is not an OH MM but is similar to Arp220 
in its properties. Both have dominant (circum-)nuclear 
starburst activity and NGC 6240 also has a shock-heated 



molecu lar emission region l ocated between it s two AGN 
nucle i ( Tacconi et al. 19991 : Baan et al. 2007 : Iono et al. 
20071 ). The emission of the ring of starformation in 



NGC 1068 (Schinncrcr et al. 2000) also falls within the 



SEST observing beam; this would dominate the signature 
of the source and results in different ratios than for 



IRAM 30m data (jUsero et al J 12004( 1 . As a result of these 
characteristics, the three sources lie in the PDR region of 
Fig. [7K and in the XDR region of Fig. [7h. Alternatively, 
the emission from these sources could be dominated by 
low Nn emission regions, where the theoretical line ratios 
are reversed compared to the high Njj cases (Sec. (SJ). 

The evolved starburst behavior displayed in Fig. [7] 
shows that the first-order PDR-XDR division of the three 
diagrams may not be optimal. The most reliable tracer of 
PDR versus XDR behavior appears to be the HNC/HCN 
ratio, because it is least affected by late-stage starburst 
conditions. Using this ratio would also imply that most 
sources of our sample are PDR dominated starburst 
sources. 
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Fig. 10. The CN(1— 0) and CS (3—2) lines compared with other high-density tracers, a) The CN/HCN ratio versus 
the CN/HNC ratio and b) versus the CN/HCO + ratio. c) The CS/HCN ratio versus the CS/HNC ratio and d) versus 
the CS/HCO + ratio. The symbols are described with Figure 1. 



CN and CS lines - The 7V(CN)/7V(HCN) ratios 
would lie below 2.0 fo r PDRs and (far) above for XDRs 
(jMeiierink et al.l l2006h . In Figure [gib any XDR sources 
would lie (far) above the line in the two diagrams. The 
OH MM have a typical CN(l-0)/HCN ratio in the 0.4- 
0.8 range, except for Arp220 with a value of 1.7. The 
PDR-XDR division appears well situated in the diagram. 

The 7V(CS)/7V(HCN) ra tio should be > 1 at high col- 
umn densities for both XDRs and PDRs (Fig. Et). The 
observed CS(3— 2)/HCN(l— 0) ratio appears to increase 
with increasing Lfir, which suggests a correlation with 
column density. 

9.2. High-density tracer ratios versus each other 

Tracers HCN, HCO+, and HNC - The line ratios of 
high-density tracer molecules HCN, HCO + , and HNC may 
be compared to each other in order to discern collective 
behavior (Fig. [5]). There is a significant spread in data 
points and a separation of a group of mostly OH MM 
from an extended group of non— OH MM. The discrepant 
OH MM point is again IC 694. 

The diagrams of Fig. [5] may be interpreted in the 
framework of dominant PDR and XDR characteristics of 
the emission regions using only the (more dependable) 
HNC/HCN ratio (see Section [5]). The data points have 
extended (slanted) distributions in all three frames with 



most sources lying on the PDR side of the three dividing- 
lines. 



Th e dependence 



HCO + (|Meiierink et al 



on d ensity 
20071) . In 



is 
the 



strongest for 
top-left frame, 
the diagonal structure represents a change from higher 
(bottom-right) to lower (top-left) densities with Arp220 
and other OH MM lying at the high density extreme. 
HNC depletion and HCO + enhancement due to evolved 
starburst environmental effects (see Sec. H]) may account 
for the bunching of data points, and the fact that some 
sources have moved close to the XDR-PDR boundary. 

Tracers CN and CS - The CN/HNC and CN/HCO+ 
versus CN/HCN distributions show a significant spread 
with the OH MM concentrated at lower CN/HNC ratios 
(Fig. [TUkV The first order diagnostics discussed in Sec. 8 
indicates that all sources have a CN/HCN ratio close to 
unity, which is consistent with PDR environments. There 
is a hint of a gap in the source distribution at CN/HCN « 
1. The slant in the data distribution in the bottom frame 
of Fig. [TUk (as compared with the top frame) results from 
the HCO + enhancement for part of the sample. 

The CS diagrams (Fig. [TOb) display a clear extended 
structure with the OH MM showing high ratios against 
HNC, HCO+ and HCN. Since CS is to first order a col- 
umn density tracer (see Sec. [5] and [5]), this linear variation 
in the relative strength of CS is a consequence of seeing 
more of the highest density components with a higher col- 
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symbols are described with Figure 1. 



uma density. ULIRGs/OH MM with high column densities 
have high ratios and moderate (low) power FIR sources 
like NGC 1808, M83, and Maffei2 have smaller values. 

Comparisons with CO - A diagnostic diagram 
may be constructed using the ratios of tracers with 
CO (1—0) (Fig. [TT]). These diagrams also display a sig- 
nificant (diagonal) spread in data points and a separa- 
tion between OH MM and the rest of the data points. 
The HCO + enhancement and the less apparent HNC de- 
pletion may account for the curved distributions of the 
bottom-left and top-right frames. As presented in previ- 
ous sections, the ratios of high-density tracers with CO 
would be evolutionary indicators that reduce in value 
during the course of the outburst. Depending on their 
LpiR(max), the sources during early evolutionary stages 
have a starting point further towards the top-right side 
of each frame, while more evolved sources or those with 
smaller LpiR(max) will appear more towards the bottom- 
left side. The relative locations of (late-stage) M 82 and 
(early-stage low L F i R (max)) NGC 253, and ULIRGs/OH 
MM Arp220 and Mrk231 are consistent with this pic- 



ture. While there is no one-to-one interpretation of rel- 
ative positions of individual data points, the collective 
paths would support FIR-related evolution from higher 
ipiR and PDR— dominated sources in early stages of evo- 
lution to XDR— like sources at lower Lfir (see Figs. |H1 
[TOl ITTl and ITS]) . Detailed modelling needs to confirm the 
evolutionary paths followed by the nuclear activity having 
different initial conditions. 



9.3. Diagnostics of the nuclear energy source 

The variation of the line intensity ratios has also been 
connected to the AGN or starburst nucleus (SBN) na- 



ture of the nuclear pow er sources (jKohno et alj 12001 



Gracia-Carpio et al. 2006f ). Both authors suggest that the 
data points with HCO+/HCN < 0.55 would indicate AGN 
activity and that higher values indicate a composite AGN 
nature. 

The (log- log) diagram of HCO + /HCN versus the evo- 
lutionary parameter of the HCN/CO ratio (Fig. [T2"|) shows 
a continuous distribution for our larger sample (see Sees. 
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and [7J) that c onfirms the distribution of data points 



prese nted earlier (jKohno et alj|2001i IGracia-Carpio et al 



2006). The central ridge that accounts for most of the 
data points shows an inverse-linear relation between the 
two ratios, which suggests a general depletion of HCN rel- 
ative to CO and HCO + during the course of the outburst. 
In our evolutionary scenario, a source moves from (bot- 
tom) right to (upper) left in the diagram during the decay 
segment of the nuclear outburst. M82 represents a later 
stage of starburst evolution and is located at the upper 
left, since its late-PDR nature also enhances HCO + . OH 
MMs and ULIRGs may start at the right part of Fig. [121 
while less dominant starbursts have a starting point more 
at the center or even on the left hand side of the diagram. 
NGC 253 is representing such a case. 

First order diagnostics (Sec. [8]) relates the 
HCO+/HCN ratio to PDR and XDR signatures rather 
than to the nature of central power source. We hereby 
assume that high Ah PDRs dominate the emission sig- 
natures and produce the low values for the HCO + /HCN 
ratio, rather than low A H XDRs. Of the four OH MM 
sources with the lowest (< 0.55 or log < -0.26) ratios 
(Arp220, IRAS 11506-3851, IRAS 15107+0724, and 
Mrk273), the radio classificatio n only indicates a p ossibl e 



AGN in IRAS 15107+0724 (|Baan fc Klocknerl 120061 ) 



NIR classification sug gests an AGN powe r contribution 



of 50% in Mrk273 (jGenzel et all Il998fh Both these 
classifications suggest a significant AGN presence in 
Mrk 231, which has a ra t io her e of 0.72 (higher than in 
e.g. IGracia-Carpio et al.1 l2006h . Observed HCO+/HCN 
values are also higher than 0.55 for the starb urst plus 
AGN source NGC 1068 f0.9: lUsero et~aD 12004^1 and the 



nearby st arburst-dominated sources NGC 253 (0.9) and 
M82 (1.6: lFuente et al.ll2005l ). 

Therefore, the predicted properties of the nuclear ISM 
during the evolution of the starburst may well account for 
the data points in Fig. [12] and particularly the (inverse- 
linear) central ridge. If the observed trend were affected 
by the presence of an AGN, its influence on the excita- 
tion environment of the nucleus would have to be further 
confirmed by more detailed modelling. The (circum-) nu- 
clear starburst would easily outshine the compact AGN 
region and dominate the emission signature in low spatial 
resolution data. 

10. Summary 

Molecular line emissions in FIR-luminous galaxies are 
tools for multi-dimensional diagnostics of the environmen- 
tal parameters in the nuclear ISM and the heating pro- 
cesses resulting from the nuclear activity. A total of 118 
line detections made with the 30-m Pico Veleta and the 
15-m SEST telescopes for a total of 37 sources have been 
complemented with partial records published for another 
80 sources. 

A proper evaluation of integrated emission lines and 
their diagnostic line ratios can be achieved after synthesis 
and modelling of a representative description of the inte- 



grated nuclear molecular medium under a wide variety of 
AGN- and starburst-related circumstances. The molecu- 
lar information in this paper has been diagnosed to first 
order using modelling results for nuclear emission regions 
presented in the literature. 

The collective behavior of line luminosities and line ra- 
tios of the low-density and high-density tracers presents a 
consistent picture of the molecular medium in the nuclear 
regions of ULIRGs. The high-density tracers represent the 
molecular medium in the regions where star formation is 
taking place, and the low-density tracer represents the rel- 
atively unperturbed larger-scale molecular environment. 
The tracer luminosities increase roughly linearly with FIR 
luminosity but they show significant scatter in data points 
due to physical processes in the nuclear region, the chem- 
ical history, and the relative age of the nuclear activity. 
The luminosity of high-density tracers varies linearly with 
Lfir except for the CS(3— 2) luminosity, where a steeper 
relation is suggestive of a different excitation dependence. 
The CO(1-0) and CO(2-l) lines for these high FIR lu- 
minosity sources have a slope less than unity and are (al- 
most) consistent with emp irical evidence o n the relation 
of SFR versus gas content ( Kennicuttl 19981 ). 

First order diagnostics of the nuclear ISM can be 
based on the collective behavior of high-density tracers 
HCN, HNC and HCO + showing significant differentiation 
and systematic changes. The line strengths of these three 
molecules relative to the CO (1—0) line show a significant 
range of 0.03 to 0.23 with a distinct dependence on Lfir- 
This dependence has been interpreted in terms of an evo- 
lutionary model where the depletion of molecular gas de- 
pends on the consumption and destruction of the high- 
density gas by the ongoing star formation process. Only 
partial diagnostics could be done using the available CN 
and CS data, which complements the diagnostics of the 
three other high-density tracers. 

The emission line ratios of the three high-density trac- 
ers of the galaxies show a structured distribution that 
fills selected parts of the parameter space. Interpreting 
these distributions has been done using modelling col- 
umn densities that are characteristic of PDR- and XDR- 
dominatcd nuclear environments. OH MM and other 
powerful ULIRGs are mostly characterized by PDR- 
dominancc. The other end of the data distribution is char- 
acterized by mixed PDR-XDR and XDR-dominated envi- 
ronments. (U)LIRGs represent the phases of nuclear evo- 
lution that generate the highest FIR luminosities and most 
rapidly deplete the dense molecular component. Using the 
HCN/CO ratio as an evolutionary indicator, the distribu- 
tions of data points may reflect evolution from PDR- like to 
more XDR-like nuclear ISM properties during the course 
of the outburst. 

The HCO+/HCN and HNC/HCN ratios serve as in- 
dicators of environments with shocks and non-standard 
heating and with HNC depletion and HCO + enhancement 
in a fraction of the sources. The HCO + /HCN ratio also 
serves as a density indicator and would be higher under 
PDR circumstances and lower in XDR environments. The 
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observed trends in the HCO + /HCN ratio could result from 
the (indirect) influence of an AGN in the nucleus. More 
likely the trends in the relative abundance of HCN as com- 
pared with other constituents result from evolution of the 
nuclear environment during the SBN activity. 

The detailed interpretations of multi-line multi- 
molecule emission line behavior and their relation with 
the local environment require detailed modelling of the 
physical parameters of the environment together with the 
excitation, the chemistry, and the radiative transfer for the 
molecular constituents. Hereby the integration of higher 
level transitions of the molecules is needed to determine 
specific excitation temperatures and densities. A compar- 
ison with the properties of Galactic emission regions will 
emphasize the effect of scale size on the line ratios for 
tracer molecules. Further studies are underway to connect 
the FIR signature a nd global heating scenario of ULIRGs 
( Loenen et ahll2006h with modelling of the molecular emis- 
sions under varying conditions in a nuclear starburst. The 
emission scenarios for XDRs and PDRs establish the con- 
nection between star formation and other sources of ex- 
citation and the integrated emission line parameters ob- 
served in extragalactic sources. 
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Fig. 2. a) Integrated HCN(l-O), b) integrated 
HNC(l-O), and c) integrated HCO+ (1-0) line lu- 
minosity versus FIR luminosity for the total sample. 
There are 84 reliable data points for HCN, 27 for HNC, 
and 42 for HCO + . The slope of the fitted line (s) and 
the correlation coefficient (r) are shown in the lower right 
corner of each frame. The displayed upper limits have 
not been used in the fits. The symbols are explained with 
Figure 1. 



Fig. 3. a) Integrated CN N = (1-0), b) integrated CN 
N = (2—1), and c) integrated CS (3— 2) line luminosity 
versus FIR luminosity for the total sample. There are 19 
reliable data points for CN(l-O), 13 for CN(2-1), and 
17 for CS(3— 2). The slope of the fitted curves (s) and 
the correlation coefficient (r) are shown in the lower right 
corner of each frame. No upper limits have been used in 
the fits. The symbols are described with Figure 1. 
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Fig. 4. The fraction of the high-density com- 
ponent relative to the low-density component, 
a) Integrated HCN (1-0)/CO (1-0), b) inte- 
grated HNC(1-0)/CO(1-0), and c) integrated 
HCO+(1-0)/CO(1-0) line ratios as a function of 
FIR luminosity. The symbols are described with Figure 1. 



Fig. 5. The fraction of the high-density component 
relative to the low-density component, a) Integrated 
CN(1-0)/CO(1-0), b) integrated CN (2-l)/CO (1-0), 
and c) integrated CS (3— 2)/CO (1— 0) line ratios versus 
FIR luminosity. The symbols are described with Figure 1 . 
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Fig. 6. Simulations of characteristic behavior of the high- 
density component during an FIR outburst. Labels (a) 
and (b) indicate the rapidly rising and the slower decay- 
ing segments of the outburst, respectively. Upper diagram: 
Three Lfir light curves versus time on a log-linear scale. 
The time is in units of T as used in Eq. 1. Middle diagram: 
The associated luminosity of the high-density (HD) com- 
ponent Lhd for the three light curves on a log-log scale. 
The parameters of the curves as described in Sec.[7]are: 1) 
the outer curve with initial (3q = 0.18, log LpiR(max) = 
12.5, and a depletion rate of T — 0.9 during the course of 
the outburst; 2) the middle curve 0.135, 12.0, and 0.5; and 
3) the inner curve 0.09, 11.0, and 0.4. Lower diagram: The 
time-evolution of the ratio of high- and low-density (LD) 
components for the three FIR light curves of the above 
frames. 
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Fig. 7. High-density tracer line ratios versus FIR lumi- 
nosity, a) the integrated HNC (l-0)/HCN (1-0) line 
intensity ratio, b) the integrated HCO + (1-0)/HCN (1-0) 
line intensity ratio, and c) the integrated 
HNC(1-0)/HCO + (1-0) line intensity ratio. The dotted 
lines mark the predicted range of XDR-PDR excitation 
using first order modelling diagnostics. In the upper 
frame the PDRs lie below the line and the XDRs above. 
These labels have been omitted in the two lower frames 
(where X would be in the upper range of frame b) and 
in lower range of frame c)) because non-standard heating 
effects on HCO + and HNC (see section l97l|) cause a shift 
of the data points and makes labelling non-appropriate. 
The symbols are described with Figure 1. 
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Fig. 8. CN/HCN and CS/HCN line ratios versus FIR 
luminosity, a) The integrated CN (l-0)/HCN (1-0) line 
intensity ratio, b) the integrated CN (2-l)/HCN (1-0) 
line intensity ratio, and c) the integrated 
CS(3-2)/HCN(l-0) line intensity ratio. The sym- 
bols X and P and the dotted lines mark the range of 
values for XDRs and PDRs. For the CS/HCN diagram 
the dotted line divides the range for low- and high-column 
densities (L & H). The symbols are described with Figure 
1. 



